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Within two Higgs doublet models, it is possible that the current vacuum is not the global min¬ 
imum, in which case it could possibly decay at a later stage. We discuss the tree-level conditions 
which must be obeyed by the most general scalar potential in order to preclude that possibility. 
We propose a new procedure which is not only more general but also easier to implement than the 
previously published one, including CP conserving as well as CP violating scalar sectors. We illus¬ 
trate these conditions within the context of the Zo model, softly broken by a complex, CP violating 
parameter. 

PACS numbers: 12.60.Fr, 14.80.Ec, 14.80.-j 


I. INTRODUCTION 

The discovery of a fundamental scalar particle mm prompts the search for other fundamental scalars. The simplest 
theoretical possibility consists in adding a second doublet to the Standard Model (SM), forming a two Higgs doublet 
model (2HDM) - for reviews see, for example, Refs. [E, 4j. 

There are several theoretical constraints that a viable 2HDM scalar sector must obey: it must be bounded from 
below [5], and it must conform to perturbative unitarity IBHS1. In addition, the scalar potential of the 2HDM may 
have simultaneously two neutral minima, both CP conserving or both CP-violating mm- In those cases, if one 
were in the local metastable minimum, there would be the possibility of decaying later into the global minimum: the 
true vacuum. We dub this possibility the panic vacuum. 

When building a model with a scalar sector allowing for several distinct minima, one is free to take a metastable 
minimum as physically acceptable, provided it is sufficiently long-lived, or to disregard it altogether. But one defi¬ 
nitely cannot neglect this issue. Typically, this would be investigated by reconstructing the potential from the input 
parameters (vacuum expectation values, masses, etc...), finding the second vacuum numerically and comparing the 
potential depths in the two vacua. This is time consuming, especially if one performs extensive scans over parameter 
space, and it could give a false positive. Luckily, within the 2HDM there is an efficient way to study these issues with¬ 
out trying to solve the minimization problem. Following the works I12H14] . a discriminant for panic vacua has been 
studied in the context of the 2HDM with a softly broken U( 1) symmetry [IS], which we denote by “17(1)-2HDM”, and 
in the context of a ^-symmetric 2HDM, softly broken by a real parameter mi which we denote by “real 2HDM”. 
The discriminant in this case is a compact combination of the parameters of the potential and the vacuum expectation 
values (vevs) calculated at a single minimum, whose positive sign guarantees tree-level stability and negative sign 
signals the presence of a deeper minimum. 

A technique applicable in the case of a scalar potential without CP violation was also presented in Ref. [IS]. The 
technique had five steps: i) perform a numerical search for the eigenvalues of a certain 4x4 matrix (A e)\ h) perform 
a numerical calculation of its eigenvectors; iii) combine the eigenvectors into a rotation matrix, used to transform the 
parameters of the potential; iv) determine whether a certain rotated quantity M 0 were positive or negative; v) apply 
the discriminant, distinguishing between the metastable vacuum and the global minimum (at tree level). 

Those studies are not applicable, for example, to the 2HDM with a complex soft breaking of the Z 2 symmetry, 
usually known as the C2HDM [201427] . There has been renewed interest in this model, in part because it is the simplest 
model allowing for the 125 GeV Higgs to have a mixture of scalar and pseudoscalar components p8H55] . Remarkably, 
all the data available is still consistent with the possibility that its couplings to the up quarks are mainly scalar while 
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its couplings to the down quarks are mainly pseudoscalar |32j . So far, no simple criteria to avoid metastable vacua 
has been devised for potentials with CP violation, either spontaneous or explicit. 

In this article, we solve that problem, thus completing the study of the 2HDM. We introduce a new discriminant 
which, using very simple criteria, distinguishes between a metastable vacuum and the global minimum. Besides 
devising a discriminant applicable to the most general 2HDM, we also improve upon previous analysis of the C7(l)- 
2HDM and the real 2HDM by presenting the discriminant in terms of physical masses, mixing angles, and couplings 
rather than the parameters of the potential. In the most general 2HDM, our procedure improves over the one suggested 
in m for the CP conserving 2HDM in two ways. First, we build it only from eigenvalues of and do not ask for 
the calculation of the eigenvectors nor for explicit basis transformation. Second, we prove that our method is safe 
even in pathological cases of potentials unbounded from below with matrix A e having complex eigenvalues. 

The paper is organized as follows. In section [IT] we present our notation, used in section [ITT] to develop the discrimi¬ 
nant D and the method which disentangles the metastable vacuum from the true vacuum. In section [TV} we discuss in 
detail the C/(1)-2HDM, the real 2HDM, and the C2HDM, presenting the discriminant in terms of physical parameters, 
such as masses and mixing angles, instead of parameters of the potential. In these cases, there is a simple condition 
forcing the potential to be bounded from below, and the method to avoid panic vacua involves a single application of 
a discriminant D. We draw our conclusions in section [v] and relegate some detailed proofs to the appendices. 


II. NOTATION 

The Higgs potential of the most general 2HDM is usually written as [1] 

V H = m 2 n \<t>i\ 2 + mlM 2 - mf 2 ^>1^2 + (mf 2 )* 4<^i 

+ “^“l0l | 4 + I 02 | 4 + A 3 |</>i|"|(/> 2| 2 + A 4 (</>{</> 2 ) ( 020 l) 
+ “^“(0102) 2 + -y-(020l)“ + 1011“ A6 0102 + Ag 020! 


+ 102 


A7 0i02 + Ay <j)\4> 1 


(1) 


where hermiticity forces all couplings to be real, except mf 2 , A 5 , Ag, and A 7 . If all couplings are real, then there is 
explicit CP conservation. 

For the study of minima, an alternative formulation in terms of quantities bilinear in the fields is useful HM3: 


V H = ~M^ + -A^rV", 


where 1 


and 


= (r 0 ,ri,r 2 ,r 3 ) = (|0i | 2 + \cj) 2 \ 2 , 2Re [</>\ 


21 m 


\<h) , 1011 2 - 1021 2 ) , 


M/j. = I — 


ni 


l 22 


, Re (mj 2 ), -Im (m 2 2 ) , -- 


m il m ‘ 


22 


( 2 ) 


(3) 


(4) 


while A^j, is a symmetric matrix. For our purposes, it is more useful to consider the mixed symmetry tensor A'',, = 
A a!/ , viewed as a matrix in Euclidean space: 


A E = A'C = - 


( §(Ai + A 2 ) + A 3 Re(Ag + Ay) —Im(A 6 + Ay) 

—Re (Ag + A 7 ) — A 4 — Re(A, 5 ) Im(A 5 ) 

Im (Ag + A 7 ) Inr(A 5 ) -A 4 + Re(A 5 ) 


|(Ai — A2) 
—Re (Ag — A 7 ) 
Im (Ag — A 7 ) 


\ 


V “ 5 (Ai~A 2 ) 


(5) 


—Re (Ag — A 7 ) Im(Ag —A 7 ) — |(Ai+A 2 ) + A 3 / 


1 We adhere here to the original notation of |12H14I inspired by relativity, where it is the contravariant vectors which define the spatial 
components with the positive sign: = (Vfi, Vi, V 2 , V 3 ). This is the opposite of the notation used for contravariant vectors in Ref. m■ 
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Of course, one can change the basis of the fields by 


# = U ij( t> 3 . ( 6 ) 

Henceforth, using a roman lowercase letter in a four dimensional vector refers implicitly to the spatial components: 
rjfe, k = 1,2,3. Choosing the matrix U in U( 2) guarantees that the kinetic terms retain their canonical form. In 
that case, ro is invariant and suffers a 0(3) change. However, as hinted in [TDJ QTj and explored extensively in 
Refs. [MTU, if one wishes to study exclusively the properties of the potential, one can advantageously take any 
general transformation U (unitary or not), in which case r M suffers a .50(1,3) change. The advantage of these more 
general transformations is that, for potentials bounded from below (BFB) in the strong sense |T0], one can diagonalize 
A e into 


A e —> diag (Aq, Ai, A 2 , A 3 ), 


and the eigenvalues satisfy [T2] 


A 0 >0, A 0 > A fe . 


(7) 

( 8 ) 


Note that the signs of A^ are arbitrary. 

In the light of future discussions, it is instructive to see what changes within this formalism if the potential is 
unbounded from below. Although this situation is unphysical, one might accidentally run into it when scanning over 
the parameter space. In order to be sure that the discriminant we develop below for the general 2HDM does not 
produce false positive results, we must extend the bilinear formalism to this situation. 

First, since A^ is not symmetric, it might happen that some of its eigenvalues are complex. There can be only 
one pair of complex and mutually conjugate eigenvalues, see Appendix [C] This leads to a potential unbounded from 
below, as proved in the Erratum of Ref. m- Note that in this case we can still have a local minimum; the simplest 
example being the inert 2HDM with A 2 < 0. Second, if all eigenvalues are real and some of them are degenerate, it is 
possible that A E is not diagonalizable at all. This can happen for a potential whose quartic part has flat directions, so 
that it is not bounded from below in the strong sense. Third, even if A E is diagonalizable 0 , it might happen that its 
eigenvalues do not satisfy (| 8 |, in which case the potential is also unbounded from below; yet a local minimum can exist. 
Thus, a discriminating procedure must be guaranteed to cut off these cases without referring to the diagonalization 
of A e . 

After spontaneous symmetry breaking (SSB), the fields acquire vacuum expectation values: 

^ = 71 ( v 1 ) ’ ^ = ^>(v 2 e iS ) ’ (9) 

where, without loss of generality, we have taken (</>i) to be real. Thus, 

r M = (r 0 ,ri,r 2 ,r 3 ) = (r^) = \ ( v 2 , 2 ;qu 2 cos <5, 2uiu 2 sin<5, v\ - u|) , ( 10 ) 

where v = \/v\ + v\ = (v^G ^) -1 / 2 = 246 GeV. Then, v\ and v 2 depend only on tan/3 = v 2 /v\. Notice that refers 
to bilinears in fields, while r^ refers to their vevs. Given our definitions, the allowed satisfy 

ro > 0 , > 0 , ( 11 ) 

corresponding to the forward lightcone. The SU( 2) xU( 1) symmetric vacuum lies at the apex, the surface corresponds 
to neutral vacua, while any point in the interior of the lightcone represents charge breaking vacua. 


III. DISCRIMINATING THE GLOBAL MINIMUM FROM METASTABLE VACUA 

Phenomenological analyses of 2HDM usually start with the following procedure. One assigns values to physical 
observables, such as masses, mixing angles and couplings, which obey current (and simple) experimental constraints. 
For example, one imposes that one of the Higgs masses equals 125 GeV. The fact that masses squared are (of course) 
chosen positive, implies that that point in parameter space corresponds to a minimum; it may be local or global, but it 
is guaranteed to be a minimum. Further, since either or both vevs are nonzero, one is guaranteed to be in a minimum 
which breaks 517(2) x 17(1). But one is still unprotected against unpleasant situations where one is sitting in the 
panic vacuum or when the potential is unbounded from below. From the input values, one extracts the parameters 
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of the potential. One can then impose theoretical constraints: boundedness from below, perturbative unitarity, and 
avoidance of the panic vacuum. 

When A 6 = A 7 = 0 , A E is block diagonal, its eigenvalues can be found analytically, and the conditions of bounded 
from below in the strong sense can be written simply in terms of the parameters of the potential as: 

Ai > 0 , A2 > 0 , > —A3, \J A1A2 > |A5 1 — A3 — A4. ( 12 ) 

When either Aq or A7 differ from zero [12 , one needs to construct the matrix Ag, diagonalize it numerically - c.f. 
Eqs. 0 and 0 - and then impose Eqs. | 8 ]). It turns out that, for a general potential, this part of the procedure 
is more difficult than it seems. The point is that, even after the diagonalization leading to four eigenvalues A a 
(a = 0,1, 2, 3), one must identify which one of the four is Ao. Ref pj2 addressed this issue in the context of the CP 
conserving two Higgs doublet model by resorting to the identification of eigenvectors, a rotation matrix, and some 
rotated vectors. 

Here we propose a procedure valid for the most general potential and involving only the eigenvalues. Imagine that 
we have identified the four eigenvalues A Q , that they are real, and we wish to know which one to ascribe the subindex 
0. We construct the projection operators: 


pa 


n 

0 


1 

Aq. a^ 


(A e — Apt). 


(13) 


Then, 


Sol 


sign 



(14) 


will be positive for only one value of a; the corresponding eigenvalue is the time-like A 0 . This assertion is proved in 
Appendix [A] 

The method proposed here to avoid panic vacua uses exclusively the stationarity conditions and the diagonalization 
of A e already needed to impose boundness from below. The former can be obtained by minimizing the auxiliary 
potential 


V = V-\ Cr^r M , 

with respect to and C, yielding 

A^ v r v — AP = C r M . 


(15) 


(16) 


Any component of this equation can be used to determine (. Typically, equating ( obtained from two different 
components yields some mf - (in M M ) in terms of the quartic coefficients A*, (in A^ v ). The quantity ( calculated at 
any neutral stationary point has a direct physical interpretation mm- the charged Higgs mass squared is equal to 
£v 2 . Our new discriminant is given by 


D = -det (A E - Cl), (17) 

where 1 is the four dimensional identity matrix. Writing A E in the diagonal basis, 

D=( Ao-C)(C-Ai)(C-A 2 )(C-A 3 ). (18) 

With these definitions, we introduce the following Method: 

1. Determine C from Eq. ©. 

2. Determine D from Eq. @- 

(a) If D > 0, then we are in the global minimum. 

(b) If ZA < 0, then we must continue. 

3. Find the eigenvalues of A E - If some of them are complex, discard the point. If all of them are real, use s a in 
Eq. (14) to identify A 0 . 

(a) If A 0 < A fc for some k = 1,2,3, then the potential is not bounded from below and the point must be 
discarded. 
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(b) If Ao > Afc for all k = 1,2, 3, then we must continue. 

i. If £ > A 0 then we are in a global minimum (there is no other minimum). 

ii. If £ < Ao then we are in the panic vacuum (metastable minimum). 

As mentioned, in cases other than softly broken Z 2 the analysis of boundedness from below requires the diagonalization 
of A e, while D in Eq. © can be calculated in any basis. And, in many cases of interest it will turn out that our 
point has D > 0. This explains why we compute D before checking whether the potential is bounded from below. 
When D > 0, we are in the global minimum and we automatically know that the potential is bounded from below. 
There exists no other situation in which D > 0 is compatible with a minimum. This is proved in Appendices [B] and 
0 Thus, if D > 0. metastability is avoided and the phenomenological analysis of this point in parameter space can 
continue immediately. The diagonalization of A# is only performed when D < 0. This procedure solves completely 
the identification of the global minimum in the most general 2HDM. 

Since ( = m 2 H ±/v 2 , we can also make another curious observation. Suppose we know completely the quartic part of 
the potential and we know that it is bounded from below. Then, the local vs. global minimum ambiguity is decided 
just by the value of the charged Higgs mass through Eq. ©■ 


IV. APPLICATION TO A SOFLTY BROKEN Z 2 POTENTIAL 

This section is dedicated to potentials with a softly broken Z 2 symmetry, where Xq = A 7 = 0 in Eq. 0 - We allow 
m \ 2 to be complex (C2HDM) or real (real 2HDM). In addition, we may have A 5 = 0 (t/(l)-2HDM). In such cases, 
A^; in Eq. |5]) is block diagonal and its eigenvalues are easy to find: 

Aq = ^ (^3 + V^) , A 3 = - (a 3 — , 


A-l — — ^ + 1^5 I) ' 


^2 — “ 1^51 


This formulas lead from Eqs. ([ 8 ]) to Eqs. ( |T2| . Moreover, using a basis where both vevs are real 2 , we find 

A 4 + Re(As) Re(?n 2 2 ) 


c = - 


2 

A4 + A5 
2 

l H± 


ViV 2 


+ 


m 


12 


V\V 2 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 


In going from the first to the second line we have used the the stationarity conditions. Some relevant additional 
formulas will be presented in section [IV A[ 

Since both Eqs. 0 and (jTsj) depend on A a , there is an intimate connection between boundedness from below and 
the discriminant for panic vacua. In cases other than soflty broken Z 2 (real or complex), including in particular the 
most general CP conserving potential, there is no simplified compact form to assure that the potential is bounded 
from below written in terms of the original parameters of the potential A 1 -A 7 . Thus, as proposed in our Method 
of section III it is best to start by computing D (a determinant which can be calculated in any basis) and only 
later perform the diagonalization of A^ which enables the application of Eqs. ([ 8 ]), guaranteeing that the potential is 
bounded from below. In contrast, for soflty broken Z 2 potentials (A 6 = A 7 = 0, mf 2 real or complex), the bounded 


from below conditions can be applied directly in the form of Eqs. (12). In that case, one can start by imposing these 
conditions and later study the discriminant 


D = 


D 


Ao - C 


= (C-a 1 )(C-a 2 )(c-a 3 ) 


Y< ± . A 4 y 

|A 5 | 2 ' 


Iv v 2 + 2 ) 

4 



Then: 


TO ff± VA 1 A 2 — A 3 

Q 


( 23 ) 


2 One can always find such a basis. The phase will appear as part of m 2 2 and A5. 
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We are in a global minimum if and only if D > 0. 


We stress that the use of D is only relevant after one has imposed boundedness from below. This is simple for softly 
broken Z 2 models, but requires the diagonalization of A# otherwise. In those more general cases, one would be better 
served using D before diagonalization of A#; if D > 0 the point corresponds to the global minimum (the potential 
is guaranteed to be bounded from below) and the phenomenological analysis may continue; only if D < 0 should the 
diagonalization of proceed in order to determine whether the point corresponds to a global minimum or not, along 
the lines of the Method of section IIIII 


Using Eqs. (191 and (20 1 , it is easy to write both D and D in terms of the parameters of the potential. We will 
shortly use this to describe D in terms of physical parameters. However, the presence of the square root vWlX? in 
Ao makes it impossible to write simple expressions for D in terms of physical parameters. In contrast, it is possible 
to write D in terms of physical parameters. This is what we turn to next. 


A. The C2HDM 


We consider here the C2HDM [2011281 . for which X e = A 7 = 0 and arg(As) 7 ^ 2arg(m? 2 ), ensuring that there is 
explicit CP violation in the scalar potential. We follow the notation of Refs. E3E2- In particular, a basis is chosen 
such that both vevs are real - S = 0 in Eq. ([ 9 ]). The minimization conditions are: 

— 2 to?! = -2 Re (to? 2 ) — + Ai v 2 + A345 v\ 

i'i 

—2 to? 2 = —2 Re (to? 2 ) — + A 2 v\ + A345 v\ 

V2 

2 Im (m? 2 ) = t)i«2 Im (A5), ( 24 ) 

where A345 = A3 + A4 + Re (A5). We may parametrize the original fields as 

cpG+ - sgH+ \ / s p G+ + c p H+ 

+ rn + icpG 0 - ispr] 3 ) J ’ ~ \ ^{v 2 + rn + ispG 0 + icgr] 3 ) 

where G + and G° are the Goldstone bosons and H + is charged Higgs with mass rn H ±. Thus forth c ang ie (s an gie) refers 
to the cosine (sine) of that angle. The neutral mass matrix Ai 2 is diagonalized by the orthogonal transformation 




( Vi\ 

= R I m I , (26) 

such that 

R A4 2 R t = diag (to?, to?, to?) , (27) 

where the neutral scalar masses are ordered as toi < m 2 < m 3 . We parametrize 

( C1C2 S1C2 S2 

-(C1S2S3 + S1C3) C1C3-S1S2S3 c 2 s 3 

-C1S2C3+S1S3 -(C1S3 + S1S2C3) c 2 c 3 


the matrix R by 


( 28 ) 
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where i is short for a,. The limit of the real 2HDM is obtained by setting a 2 = 03 = 0 and a = a\ — n/2. The 
parameters of the potential may be traded for physical observables through 


v 2 Ai = — 


1 


v 2 A 2 = - 


cos 2 /3 
1 

sin 2 /3 


v A3 


A 3 — 


[-ml c 2 c 2 - m|(c 3 Si + Cis 2 s 3 ) 2 - m 2 (cic 3 s 2 - S1S3) 2 + p 2 sin 2 /3] , 
[-m 2 S 2 c2 - To 2 (cic 3 - SiS 2 S 3 ) 2 - TO 2 (c 3 SiS 2 + CiS 3 ) 2 + / I 2 COS 2 /3] , 

Sm/3cOS/3 I-' 1 ° 2 1712 ( S 2 S 3 c 3l + m 3 ( S 2 C 3 s 3)) c l s l 

+ (to 2 - to 2 )(c 2 - s 2 )s 2 c 3 s 3 ] - X 2 + 2?nX± , 

2 \ 2 2,/ 2 2 1 2 2\ 2 1 2 q 2 

i? A 4 = s 2 + (rn 2 s 3 + m 3 c 3 )c 2 + /i — Zm H ±, 

2t>/\\ 22 / 2 2 1 2 2\ 2 1 2 

v Re(A 5 ) = -m 1 s 2 -(m 2 s 3 + m 3 c 3 )c 2 +p , 

2 

u 2 Im(A 5 ) = -c 2 [(-to 2 + to 2 s 2 +to!c 2 )cis 2 + (to 2 - to 2 )sis 3 c 3 1 , 

sm p 


with 


k = 


V\ V2 


Re(TO? 2 ). 


The observables are toi, to 2 , to#±, 01 , 2 , 3 , /3, and Re(m 2 2 ), while to 3 is determined from 


TOo = 


to 2 i?i 3 (i?i 2 tan/3 - i? n ) + m| R 23 (R 22 tan/? - i? 21 ) 
-^ 33(^31 — -R 32 tan/3) 


(29) 


(30) 


(31) 


We can now write the discriminant D in terms of physical parameters. Substituting Eqs. (291 into Eqs. (19) and 

1 


(20), and those into Eq. (18), we find 

D = 


8v 8 si c 2 


(-ai p 2 + 61 ) (a 2 p 2 - 2 b 2 ) , 


(32) 


where 


2 r_2 2 1 / 2 2 1 _2 21 21 

Ol = Sp [TO!S 2 + (TO 2 S 3 +TO. 3 C 3 ) C 2 J , 

61 = c 2 [cis 2 (-to 2 + TO 2 Sg + m\c\) + S 1 S 3 C 3 (to 2 - TOg)] ' , 
a 2 = 2 to 2 c 2 c 2 1+/3 + (to 2 + to 2 ) (1 - c^c 2 1+/3 ) 

+ (to 2 - m\) [cos (2a 3 ) (s 2 i+/ 3 - c 2 1 +/3 s 2 ) + sin (2a 3 )s 2 sin (2c*i + 2/3)] , 
o 2 = (to 2 c 3 + m 3 s 3 ) TO!C 2 + to 2 to 3 s 2 . (33) 


B. The real 2HDM and the £/(l)-2HDM 


We now consider a real softly broken Z 2 potential, where A 6 = A 7 = 0, and A 5 , to 2 2 , and the vevs are real. In this 
case, CP is conserved, 773 = A is CP odd, while the neutral CP even scalars are denoted by H (heavy) and h (125 
GeV). In terms of the notation in the C2HDM, the angles become a 2 = a 3 = 0 and a = oq — n/2. In this case, the 
expression for D in terms of physical parameters simplifies considerably and can be programmed directly into any 
phenomenological analysis. We find 


4 ?; 8 r 3 „3 
4 v c p s p 


D = to 2 2 



™ 2 12 

m h m H S pCp 


(' m H c a+p 


*1 2 \ 
™>h*OLM 3) 


(34) 


The discriminant introduced in reference [T9] for the case of the real 2HDM is 

Dz 2 = m i 2 ( m ii — /c 2 TO 22 )(tan/3 — k), 


( 35 ) 

















where 


k = 



(36) 


This quantity can also not be written directly in terms of physical parameters. The problem, as explained above, is 
related to the presence of \A 1 A 2 in the A 0 of Eq. (19), which is the origin of k in Eq. (35). 

Although this expression is written in terms of the parameters of the potential and not the physical observables, it 
can be recovered in our formalism in the following way. We first calculate D using (19), 


m i 2 fc4 (ra| 2 ) 2 ~ ( m n) 2 
V\V 2 v 2 uf — k 4 vf 


(37) 


and then relate it to Dz 2 via 


D = X 


2 Vl ^12 - A5U1U2 

2 uf (Ai«f — A 2 u|) 2 


(—mh — fc 2 TO2 2 )(tan/3 + A:)(tan 2 j3 + k 2 )Dz 2 , 


(38) 


where each factor is positive definite. Combining the pre-factors in Eqs. (35) and (38), we that D depends only on 
/c 4 and, thus, we understand why it can be written in terms of physical parameters in a simpler way. 

Things get even simpler in the case of the potential with a softly broken U( 1) symmetry. The discriminant for this 
case was introduced in Ref. HI] as 


Du( 1 ) = — fc 2 m 22 )(tan/3 — k). 


(39) 


For this model, 2 m \ 2 = m^s 2 /3 > 0, and the discriminant may be taken as 




D = 1 - 


m. 


b H 


L H 


(11 , 11 \ 
\ m H C a+P + m h S a+p)- 


(40) 


V. CONCLUSIONS 


When studying models with two Higgs doublets, it is possible that the vacuum chosen is not the global one, and that 
there is another vacuum lying below. In that situation, the vacuum is metastable, there could be a later transition 
into the global vacuum, and we dub this situation the panic vacuum. We have developed a method to avoid panic 
vacua, involving the discriminant D in Eqs. ~(p~8|), which is applicable to any 2HDM potential. In particular, we 
proved that D > 0 guarantees that we are staying in the global minimum and, automatically, the potential is bounded 
from below. If D < 0, more steps are needed to discriminate between the global minimum from a panic vacuum or 
a potential unbounded from below. Our method is not computer time consuming and it be easily implemented in 
phenomenological studies in which extensive scans over 2HDM parameter space. 

We have shown how D can be written in terms of physical parameters for the C2HDM, the real 2HDM and the 
t/(l)-2HDM, and we have shown how it is related with the discriminants previously presented for the later two cases 

QSHH. 

If the potential has a softly broken Z 2 symmetry, a simpler strategy is possible. Indeed, in that case, the bounded 
from below conditions are easy to implement - c.f. Eq. (12). Then one need only apply D in Eq. (23) and one has a 
global minimum if and only if D > 0. 

This article completes the identification of panic vacuua for all 2HDM, including the previous unsolved cases where 
CP violation is present in the scalar potential. 
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Appendix A: Proof that so > 0 


The operators P a introduced in Eq. (13), one for each eigenvalue A a of A E , obey 


P a pP = 6 a pP a , ]TP“=1, 


(Al) 


and P a is the projection operator into the subspace generated by the eigenvector corresponding to A Q . Since our 
method uses projectors after we have discarded complex eigenvalue situation, we know that all A a are real. 

The explanation of the proof is easiest in Minkovski notation, where 

= n (V - A ^), (A2) 

which are symmetric in the indices fi, v. Take now some four-vector (for example = (1,0,0,0)). The new vector 

C = {P°)u^ (A3) 

lies along the eigenvector of A 0 , and, thus, it is timelike: 

o < = ^'(P°)^(P 0 )^ = r'(p°) M> r. 


(A4) 
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Similarly, 

t = (P k ) v ^ (A5) 

for k = 1, 2, 3 lies along the eigenvector of A*, and, thus, it is spacelike: 

0 > (£ fc )^ = & (. P k (A6) 

This is true for any vector F, apart, of course, from the case when the vector is accidentally chosen to be orthogonal 
to some eigenvector. Choosing the simplest case i M = (1,0, 0,0), means that 


Sa. 


sign 



(A7) 


is positive if and only if a = 0. This completes our proof. 


Appendix B: Locating stationary points in parameter space 


Imagine that A 0 , A*,, and ( have been identified. Some ordering of A a has been found. If A 0 is the largest, we 
denote the spacelike eigenvalues by A a < A b < A c . In this appendix we show how one can determine how many 
extrema exist, depending on where ( sits with respect to the A a . We follow the analysis of Refs. nmmnsi, also 


used in m in the particular case M 2 = 0. 

The stationarity conditions (16) may be written in the basis where A e is diagonal as 


(A 0 - C)r 0 = M 0 , (A fc - C )h = M k , 


(Bl) 


where the hat in r and M emphasize the fact that these vectors are written in the basis where A^ is diagonal. Let us 
introduce the variables 


r k M k A fe - £ 

n k = —, m k = - r -, a k (Q = —- - , 

ro Mo Aq - Q 


(B2) 


where n k is a variable unit vector, m k is a fixed vector, and a k ( C) are j us f numbers which depend on £. The system 
(Bl I then takes the following form: 


a k {()n k = m k . 


(B3) 


No summation is assumed here. Take some fixed ( and all possible directions of the unit vector n k . Then, the left 
hand side of the system of Eqs. (B3) defines (the surface of) an ellipsoid with semiaxes |otfc (C) I - As C changes, this 
ellipsoid grows and shrinks in a way that can be well visualized. Whenever this ellipsoid crosses the fixed point 
(mi, m2, m3), the system (B3| is satisfied, and we get an extremum. Depending on the relations among A a and £, as 


well as on the location of m k , this can happen at most six times, yielding up to six stationary points. 

Although these properties are reasonable to assume given the work of Refs. mrnm, we motivate them here for 
completeness. We start by noting that n k = 1 and Eqs. (B3) lead to 


a Q (C) 2 + a b (0 2 + a c (C) 


= 1 , 


(B4) 


which, for some fixed £ can be viewed as the surface of an ellipsoid in the (rn a . mj,, m c ) space. We will now see 
what happens to this surface as £ increases from —00 to + 00 . Recall that we have ordered the eigenvalues such that 
A a < A b < A c . In the limit that £ —> — 00 , we obtain the unit sphere. As ( increases, the sphere becomes an ellipsoid 
whose semi-axis along m a shrinks faster than the rest. When ( = A a , a a (() = 0 and the ellipsoid becomes a filled 
ellipse on the (nib, m c ) plane, flat along m a - This process is illustrated on the left panel of Fig. [lj through a projection 
on the ( m a ,rrib ) plane. Consider a fixed point on the (m a ,rrib,m c ) space. For simplicity, we illustrate this with a 
point on the m a , m b plane, shown as a black dot in Fig. [l] During the process depicted in on the left panel of Fig. [l] 
that point is crossed only once, and, thus, Eq. (B3) is satisfied only once. This means that there can be at most one 
stationary point for a £ in this region. 

Now consider that £ G [A a , A fc ], As ( moves away from A a , a a (£) increases again from zero, while a b (C) decreases. 
This process is illustrated on the right panel of Fig. |T] through a projection on the (m a , m b ) plane. When £ = Ab we 














11 




FIG. 

A a to 


1: Projection on the (m a ,m,b) plane of the ellipsoid given by Eq. (B4 ( as £ varies from 
A;, (right panel). The black dot represents a fixed point on the (m a ,rrib) plane 


—oo to A a (left panel) and from 


reach a,b(() = 0, and the ellipsoid collapses onto the (m a ,m c ) plane. This is shown as the vertical line on the right 
panel of Fig. IT] Along this process, some points may be passed twice (as happens, for example, for the illustrated by 
the black dot)) meaning that one can have at most two stationary points. 

The analysis for £ £ [A&,A C ] follows the same lines and, again, one concludes that there can be at most two 
stationary points in this region. Finally, taking £ £ [A c , Ao], the ellipsoid grows without bound and all points in the 
(m a ,m{,,m c ) space are crossed exactly once. As shown in Refs. nnmn, larger ( correspond to smaller values for the 
potential. In our notation, if M and r' M are two stationary points with the corresponding values ( > then 


\ V (r)\-\V(r')\ = -(<;-C)r^>0. 


(B5) 


This means that the largest £ for which Eq. (B3) holds will be the global minimum (provided the potential is bounded 
from below). In the case discussed here, this would occur for ( £ [A c , A 0 ]. 

Due to the positivity of r o, the above description holds for Mq > 0. If Mo < 0, then ( can vary from Ao to +oo. If 
C starts from +oo and decreases, then the ellipsoid starts from the unit sphere and monotonously grows to infinity as 
( approaches Ao- This process covers exactly once all points in the (m a ,m6,m c ) space outside the unit sphere. We 
get a single extremum in this case, which must be the global minimum. 

Finally, if it happens that A 0 is not the largest among the eigenvalues, the above construction is still valid up to 
obvious modifications and can be used to count extrema in each region. However the potential in this case is not 
bounded from below, and this picture does not allow one to spot the presence of this fact. 


Appendix C: Proof that D > 0 guarantees a global minimum 


When building a 2HDM potential, we know by construction that we are at a minimum, and we just want to know 
whether this minimum is global. In our method, we first provide a discriminant D, whose positive sign guarantees 
that we are at a global minimum. Once again we repeat that D > 0 is a sufficient but not necessary condition for the 
global minimum, and our method contains extra steps to be checked in the case of D < 0. 

Here, we prove that this is indeed a sufficient condition. We do it by showing that all other possible stationary 
points with D > 0 do not correspond to a minimum. Since our method suggests checking D even before we diagonalize 
A e, this claim must cover all possible cases, including potentials unbounded from below. 


Non-pathological case. Let us first build the proof for the case when all eigenvalues of A^ are real and this 
matrix is diagonalizable. To distinguish a minimum from a saddle point, one must consider the Hessian 


_ 1 d 2 V 

~ 2 dipadipb ’ 


(Cl) 


where for clarity the Higgs doublets are rewritten in terms of eight real fields <^ a , a = 1,...,8. One then checks 
that, apart from would-be Goldstone modes, the Hessian is positive definite in the Higgs field space. For a neutral 
stationary point, the charged and neutral fields decouple, and one focuses on H a b in the four-dimensional subspace of 
neutral Higgs modes (labeled by a, 6 = 1... 4). 
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Within the bilinear approach, the neutral Higgs mass matrix was written in the Higgs basis in mmm- Despite 
being compact, those expressions do not provide insight on how the value of £ with respect to Ao and A/- is related 
with minimum versus saddle point assignments. To gain it, we use instead the basis-invariant approach to Higgs 
masses (34] . which allows us to switch to the A M „-diagonal basis. In this basis, the Hessian in the 4D space of neutral 
inodes takes the following form: 

Hab (R ^)aa^afiiRf3b ; 

where 

/ Ao - £ 0 0 0 \ 

o C - Ai 0 0 

o o c-a 2 o 

Vo 0 0 C-Aa/ 

1 / dr a \ 

2 \ dtp a / ' 


S 


Raa 


(C2) 


(C3) 


There are a few important remarks concerning this formula. First, all indices here refer to four dimensional Euclidean 
spaces, and H a b can be calculated as a usual product of three matrices. However, there are two different spaces 
involved here. The indices a , /3 refer to the same space of bilinears from which we erased its Minkowski space metric, 
while the indices a, b refer to the space of neutral scalar modes at the extremum. These two 4D spaces are shown in 
Fig {2] Second, the signature matrix S and the matrix A e — £1) whose determinant we calculate to get D , originate 



FIG. 2: The correspondence between the 4D spaces established by the transformation R in (C3). The neutral field space 
contains a would-be Goldstone G°, which is annihilated by R. Its complement, the physical Higgs space, shown in gray, is 
mapped onto a 3D plane tangent to the forward lightcone in the r a space at the extremum point. The direction orthogonal to 
this plane points inside the cone and it cannot be represented by any neutral scalar mode. 


from the same A ^ the former is obtained by erasing the Minkowski metric; the latter is obtained by lowering 

one index. This leads to D = det S. However, as we will see below, in our discussion we will need not only det S but 
also its full signature; defined as the number of positive and negative eigenvalues. Third, since the diagonalization 
of belongs in general, to S'0(l,3), it modifies the kinetic term of the Higgs fields, and the mass matrix is not 


equal to the Hessian (C3). However such a transformation, effectively rotating and stretching the Higgs space, leaves 
invariant the signature of the Hessian (the signs of its eigenvalues) - see proposition 1 in Ref. [T2]. Therefore, the 
signature of H a b faithfully represents the signs of the masses squared of the four neutral scalar degrees of freedom. 

The link between the signatures of H a b and of S is less trivial. The transformation matrix, R , is singular: det R = 0, 
which indicates that one of the four directions in the scalar field space ip a is, in fact, a would-be Goldstone boson. 
Orthogonal to it lies the 3D physical neutral Higgs space (thenceforth referred to as the “physical Higgs space”), which 
is mapped by R onto the 3D subspace in the r a space tangent to the lightcone at the extremum point. Therefore, 
when distinguishing a minimum from a saddle point, one must pay attention not to the full signature of S but to the 
signature of its restriction onto this 3D subspace. 

This allows us to establish the following relation between the signatures of S and H a b. 


• If S has the signature (+,+,+,+), then H a b is positive definite in the physical Higgs space, and we are at a 
minimum. Also, D > 0. 
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• If 5 has the signature (+, +, +, —), up to permutations, then H a b is not positive definite in the entire r a space, 
but it might still be positive definite when restricted to the physical Higgs space. In this case D < 0, but 
deciding whether we are at a minimum or at a saddle point requires further analysis. 

• If 5 has the signature (+, +, —, —), up to permutations, then H a i, cannot be positive definite even when projected 
onto the physical Higgs space. This results from dimension counting: in the r a space, there exists a 2D subspace 
of negative S a pr a r/s, which must intersect a 3D subspace tangent to the lightcone. So, this stationary point 
cannot be a minimum but it has D > 0. The same conclusion holds for the signature 

With this classification in mind, a minimum yielding positive D can take place if and only if 5 has signature 
(+,+,+,+). This, in turn, can happen only for A 0 > £ > A c , the largest among A*,, which automatically im¬ 
plies that the potential is bounded from below. In the previous subsection, we established that there can exist only 
one stationary point in this region of £, and that this stationary point must be a minimum. In this situation, there 
cannot be any other stationary point with £ > Aq. Therefore, a minimum with D > 0 is the global minimum. 


Pathological case. Now we turn to the case of one pair of complex and mutually conjugate eigenvalues of Ae- 
Now, A^j, cannot be diagonalized with any transformation from the group 50(1,3). Nevertheless, it can be brought, 
with the aid of an 50(1, 3) transformation, to the block-diagonal form 


/ A 0 o — Aoi 0 0 \ 

I Aoi —An 0 0 I 

0 0 A 2 0 ’ 

Vo o o a 3 / 


(C4) 


where A 2 and A 3 are the two real eigenvalues of A#. It can 
expression for the Hessian (C3) is also valid with the following 


still be used to calculate D = — det(A E — £1). 
signature matrix 5: 


The 


/ A 0 o — £ A 0 i 0 0 \ 

Ad Ah+C 0 0 

0 0 c - a 2 0 

Vo 0 0 £ - A 3 / 


(C5) 


Now, we have the following sequence of arguments. Consider the discriminant D on the entire £ axis. By construction, 
it changes sign only when £ passes through an eigenvalue, which in this case happens only at £ = A 2 and £ = A 3 . 
Since D < 0 at £ —> oo, we find that the only region where D is positive is between A 2 and A 3 . The upper block in A# 
gives the following strictly positive factor to the discriminant (£ — A 0 o)(C + An) + Aqi > 0 (otherwise, the eigenvalues 
on this subspace would be real, contrary to our assumption). The corresponding upper block in 5 has an extra minus 
sign in its second row and, therefore, it contributes one ” entry to the signature of 5. Since £ is between A 2 and 
A 3 , we get a second ” from this subspace. Overall, the signature matrix 5 possesses two negative eigenvalues and 
cannot correspond to a minimum. 



